The longest excursion of fractional Brownian motion : numerical evidence of 

non-Markovian effects 
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We study, using exact numerical simulations, the statistics of the longest excursion Z ma x(i) up 
to time t for the fractional Brownian motion with Hurst exponent < H < 1. We show that in 
the large t limit, (Z max (i)) tx Qoot where Qoo = Qoo{H) depends continuously on H, and in a non 
trivial way. These results are compared with exact analytical results obtained recently for a renewal 
process with an associated persistence exponent 6 = 1 — H. This comparison shows that Qoo{H) 
carries the clear signature of non-Markovian effects for H 7^ 1/2. The pre-asymptotic behavior of 
(lmax(t)) is also discussed. 

PACS numbers: 



Introduction. In the last few years, there has been 
a growing interest in the study of anomalous dynamics 
P, Q, where by contrast with Brownian motion, long 
range temporal correlations induce non-standard dynam- 
ical behaviors. Instead of diffusive behavior, anoma- 
lous dynamics typically displays a non linear growth 
of the mean square displacement (x 2 (t)) oc t 2H , where 
H 1/2 is the Hurst exponent. Such behaviors have 
been observed in various experimental situations includ- 
ing polymer networks [H , intracellular transport two- 
dimensional rotating flows [f| or porous glasses [y]. To 
describe theoretically such situations, various stochas- 
tic processes have been proposed and studied. Among 
them, the fractional Brownian motion (fBm) , initally in- 
troduced by Mandelbrot and van Ness 0], is currently 
playing an increasing role in this area of research. For 
instance, the fBm was recently proposed to model the 
stochastic dynamics of a polymer passing through a pore 
(translocation) 

The fBm x(t) is a Gaussian stochastic process charac- 
terized by the following two-time correlations 



(x(h)x{t 2 )} = C{t u t 2 ) 
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This implies that the incremental correlation function is 
stationary, i.e. {{x(t\) — x(t 2 )) 2 ) — |ii — t 2 \ 2H ■ For 
H = 1/2, the process x(t) is just Brownian motion (BM). 
For H < 1/2 the dynamics is subdiffusive, while it is 
superdiffusive for H > 1/2. For < H < 1, fBm is 
a non-smooth process, i.e. it has an infinite density of 
zero crossings. A relevant quantity characterizing these 
zero crossings is the distribution p(r) of the time inter- 
vals between consecutive zeros. In many cases, which 
are relevant in statistical physics, this distribution has 
a power law tail p(r) oc r -1-8 , with the persistence 
exponent [l(| [H|. A remarkable result for processes, 
Gaussian or non-Gaussian, obe ying Eq. |T|), is the ex- 
act relation = 1 - H [ij, [H Qj]. Such processes (JJ) 
appear naturally in various interesting models of statisti- 
cal physics. For instance, the fBm with H = 1/4 arises as 



fBm H=0.75 



A(t) 



FIG. 1: Intervals between zero crossings (excursions) for the 
fBm in the particular case H = 0.75, which was generated 
numerically using the Levinson algorithm. The longest ex- 
cursion Z m ax(t), studied in this paper, is defined in Eq. @. 



a scaling limit of a tagged particle in a one-dimensional 
symmetric exclusion process [l5j . It also describes the 
equilibrium temporal fluctuations of the height field of a 
d-dimensional Edwards- Wilkinson interface, and in that 
case H = (1 — d/2)/2 [3]. Another example where such 
a process as in Eq. ([T]), albeit non-Gaussian, appears is 
the Mathcron-dc Marsily model of hydrodynamic flows 
in porous media. There it describes the longitudinal po- 
sition of a particle in a d + 1-dimensional layered random 
velocity field and in that case H = max(l — d/i, 1/2) [lj|. 

For Hj^ 1/2, one expects that fBm is a non-Markov 
process [7|. However the zero-crossing properties of the 
fBm which have been investigated up to now have not 
convincingly shown the signatures of these memory ef- 
fects. For instance assuming that the intervals between 
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crossings are independent and identically distributed (re- 
newal process) yields the correct behavior for the tail of 
the distribution p(r) with 9 = 1 — H [171 ], More recently 
on the basis of a numerical computation of the correla- 
tion function of the intervals between successive zeros, 
the authors of Ref. [HI claimed that the zero crossings 
properties of fBm are actually described by a renewal 
process, which contradicts our theoretical understanding 
of this process @, [HI- One goal of the present paper 
is thus to exhibit a property of the fBm which instead 
shows that temporal correlations clearly affect the zero 
crossings properties of this process. 

To this purpose, following a recent work we study 
here the statistics of the longest excursion up to time t, 
denoted l m ax(t). For a typical realization of the fBm x(t) 
with TV = N(t) zeros in the fixed time interval [0,i] (see 
Fig. [1]) let {ti,t 2 ,-- - ,tn} denote the interval lengths 
between successive zeros and A(t) denote the length (or 
age) of the last unfinished excursion. The extreme observ- 
able we focus on is the length of the longest excursion up 
to t 

lmax(t) = max(ri, r 2 , • • • , rjv, A(t)) . (2) 

We show here that the average (Z m ax(i)) is a quan- 
tity sensitive to the non-Markovian character of fBm. In 
Ref. [Hj], it was shown that (l maK (t)) can be conveniently 
computed using the exact relation 

d(l max (t))/dt = Q(t) , (3) 

where Q(t) is the probability that the last unfinished ex- 
cursion, A(t) in Fig. [TJ is the longest one 

Q(t) = Prob[/ max (t) = A(t)} . (4) 

It was then shown that for a renewal process character- 
ized by a persistence exponent 9 < 1, one has the exact 
result [li] 

lim Q(t) = 

t—>oo 

Q* = Q*(6) = / f e , (5) 

Jo l + x"e x j Q dyy °e y 

where the superscript 'R' refers to renewal process. In 
this paper, we compute numerically Q(f) ^ for fBm de- 
fined as in Eq. ([T]) for different values of < H < 1. 
We show that, in all these cases, Q(t) — > Qoo for large 
time t, as predicted in Ref. [19| for non-smooth pro- 
cesses with < 6 < 1, which is the case for fBm with 
< H < 1. We then extract precisely the asymptotic 
value Qoq = Qoa(H) : any deviation from the value 
Qoci® = 1 — H) in Eq. ([5]) can thus be identified as 
a signature of non-Markovian effects. 

Numerical method. For the purpose of numerical sim- 
ulations we need to discretize the fBm path into a set 
of Gaussian numbers correlated through Eq. (JTJ) . Gen- 
erating a sequence x — {xi, Xi, Xt} of Gaussian 
numbers with prescribed correlations ( X i X j ) = C t j is 
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FIG. 2: Q(t) as a function of t for H = 0.1, 0.5 and H = 0.9. 
The straight lines correspond to the value of Q^(l — H) for a 
renewal process given in Eq. ([5j) . This clearly illustrates that 
the fBm is not a renewal process. 

a two step procedure: i) it is first necessary to com- 
pute the matrix A, the square root of the correlation 
C . ii) Each discrete path is then given by x = A£, 
where £ = {£i, £2, ■ ■ ■ , £t} is an uncorrelated normally 
distributed set of random variables. It is easy to check 
that paths obtained from this procedure have the re- 
quired correlation matrix: 

T 

(XiXj) = A iM A k 2 ,j{^k 2 ) = A 2 l} = dj. (6) 

fel ,fc 2 = l 

Compared to standard Brownian motion, building a fBm 
is numerically cumbersome. The Brownian motion has a 
linear cost in T and is easy to simulate paths of size 
T ~ 10 6 . For fBm, the first step involves the full di- 
agonalization of matrix Cij and limitates to T ~ 1000 
the size of the path. The second step is faster and the 
matrix-vector product needs T 2 operations. A better 
performance can be obtained for fBm thanks to the sta- 
tionarity of the incremental correlation function. The 
increments Si = xi + \ — Xi are correlated according to a 
Toeplitz matrix. For Toeplitz matrices special numerical 
methods allow to build paths without going through the 
full diagonalization of C. Here we use the Levinson algo- 
rithm which is not the fastest algorithm, but is exact for 
any value of T (for a pratical implementation see [20|). 
In this paper we show the results obtained for fBm of size 
T = 10000. 

Numerical results. We now discuss our results for Q(t) 
defined in Eq. ((4]) , whic was computed by averaging over 
10 6 samples. In Fig. [2] we show a plot of Q(t) as a 
function of t for different values of H = 0.1, 0.5 and H = 
0.9. In all these cases our numerical data are consistent 
with an asymptotic behavior 

lim Q(t) - Qoo = Qoo(ff) • (7) 

t — >oo 

We also notice that this asymptotic value is approached 
from above for H = 0.1, 0.5 and from below for H = 0.9. 



3 



0.86 



©1 



0.85 




0.75 - 
0.5 - 
0.25 - 








0.25 



Renewal process 
fBm process 



0.5 
= 1 - 



0.75 



H 



FIG. 3: Plot of Q(t) as a function of t on a linear-linear 
scale for H — 0.9. The solid line indicates the fit as in Eq. 
(fl3)l with a(0.9) ~ 0.1 and 6(0.9) ~ 0.3. Inset : Plot of 
Q(t) — Qoo{H = 0.9) (same data as in the main figure) as 
a function of t in a log- log plot. The solid line corresponds 
to a(0.9)i _6 ^ H ' : this suggests a good quality of the fitting 
procedure in Eq. (JT3}. 



FIG. 4: The triangles indicate the numerical estimate of 
Qoo{d = 1 — H), extracted from the fitting procedure in Eq. 
(|13[) . For comparison, we have also plotted Q^o(0 — 1 — H) 
for a renewal process, as given by Eq. ((5j). This plot clearly 
shows that, except for H = 1/2, the fBm is not a renewal 
process. 



In this same figure [2] we also plot, with dotted lines, the 
value of for a renewal process given in Eq. © with 
9 = 1 — H . These two values Qoo and are clearly 
different as H deviates significantly from 1/2. We have 
also checked that for all these values of H, the persis- 
tence probability Po(t) ~ t~ 1+H displays a well developed 
power law behavior for t > 1000 so that a comparison 
with Q^o(0 = 1 — H) is meaningful. Therefore we can 
conclude safely that Qoo carries the signature of memory 
effects of the fBm for H ^ 1/2. 

Although the curves for Q(t) shown in Fig. [2] indicate 
an asymptotic behavior as in Eq. J7]) with a value for 
Qoo(H) different from Q^(l — H), a precise estimate of 
this asymptotic value Qoo (70 requires more effort. The 
time dependence of Q(t) is due to discretization of the 
paths and this can be understood by studying the case of 
Brownian motion (BM) . It is well known that the density 
of zero crossings of BM is infinite : this means that if the 
BM crosses zero once, it will recross zero infinitely many 
times immediately after the first crossing. Therefore a 
proper definition of the excursions requires a regulariza- 
tion procedure. A convenient way to implement it, is to 
impose that the maximal distance from the origin during 
an excursion should be bigger than xq, where xq plays 
the role of a spatial cutoff. To compute the finite time 
behavior of Q{t) we recall that the probability po(t,Xo) 
that a BM starting from xq at t = remains positive up 
to time t (persistence probability) is given by 



' '^"»li)= Erf (Ji) = Vlr +0 (^ 



j— . (8) 
Following the derivation of Ref. [19( we can compute the 



Laplace transform Q(s) of Q(t) in the limit x\ <C t [22 



Q{s) = - 

S JO 



dx 



p (x/s)e x 



Po(x/s)e x + f* dyp (y / 's)e » 



(9) 



For small s, it was shown in Ref. [19| that Q(s) ~ 
Q*(l/2)/s where Q£(l/2) = 0.626508... 0. To un- 
derstand the effects of the discretization, one needs to 
compute the first correction to this leading 1/s behav- 
ior when s — > 0. This has to be done carefully because 
a naive expansion of the persistence probability Pa(y/s) 
beyond the leading order as suggested by Eq. {5} in the 
denominator of Eq. (JS|) yields a diverging integral over 
y. Handling this singular behavior with care yields 



QOO 

where a is given by 



Q£(i/2) , ~ 



's 



(10) 



-x„l/2 



/o (x-V 2 e-* + ^FErf(Vi))' 
Going back to real time this yields finally 

Q(t) 



7 dx = 0.23970. 




^ R (l/2) +a X /^+0(xl/t) 



(11) 



(12) 



Motivated by this result for Brownian motion (fl2l) . we 
propose to describe the data for Q(t) in Fig. @ by the 
following form 



Q(t) ~ Q 00 (i7)+a(7J)i 



-b(H) 



(13) 



In particular, from Eq. (fl^j) . one expects 6(1/2) = 1/2. 
We have checked that this form (|13[) describes very well 
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our data for Q(t) for all the values of < H < 1 that we 
have studied. In the inset of Fig. [3J we show a plot of 
Q(t) — Q oc (H = 0.9), as a function of t on a log-log scale, 
while the main figure shows a plot of Q{t) as a function 
of t on a linear-linear plot. This fitting procedure (|13|) 
hence provides a reliable way to estimate the asymptotic 
value Qoo(H). In Fig. [4] we have plotted these values 
as a function of 9 = 1 — H. For comparison, we have 
also plotted the values of (9) for renewal process (0) : 
these two curves are clearly different (except for H = 1/2 
which corresponds to Brownian motion). 

Our numerical data indicate that the exponent b(H) 
exhibits a maximum for H ~ 0.5, where 6(1/2) = 1/2. 
On the other hand, one finds that b{H) — * for H — ► 
and H — ► 1 : therefore it becomes very difficult to extract 
a reliable value for 9 close to and 1. For H = 1, the 
fBm is simply a linear function of time t, x(t) — C,t where 
£ is Gaussian random variable of unit variance. It is 
thus easy to see that Q x = 1 in that case. Although it 
is very difficult to extract a reliable value of Qoo(H) for 
H > 0.95, one expects that Qoo(H) — > 1, smoothly, when 
H — ► 1. Similarly, our data suggest that Q X {H) vanishes 
smoothly as H — > 0. Finally, we notice that a(ff ) changes 
sign for H ~ 0.7 : it is positive for 7? > 0.7 and negative 
for if < 0.7. 

Conclusion. To conclude, we have presented a numer- 
ical computation of the mean longest excursion (£ max (i)) 



for the fBm with Hurst index < H < 1. We have shown 
that (Zmax(*)) ~ Qoa(H)t for large i where Qoo{H) is a 
new interesting feature of fBm. We have also demon- 
strated that this quantity is very sensitive to temporal 
correlations characterizing this process. Therefore, at 
variance with the recent claim of Ref . [181 ] , our numerical 
results clearly show that the zero crossings of fBm can 
not be described by a renewal process. We point out that 
the quantity studied here is sensitive to the full joint dis- 
tribution of the time intervals between crossings, while 
the numerical work presented in Ref. [HI only studied 
the correlation function between two such intervals. Fi- 
nally we hope that the non trivial dependence of <5oo (H) 
shown in Fig. [4] will stimulate further analytical progress 
on the study of fBm. 
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